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Abstract 

Two methods are explained to exactly solve Maxwell's equations 
where permittivity, permeability and conductivity may vary in space. 
In the constitutive relations, retardation is regarded. If the material 
properties depend but on one coordinate, general solutions are derived. 
If the properties depend on two coordinates, geometrically restricted 
solutions are obtained. Applications to graded reflectors, especially to 
dielectric mirrors, to filters, polarizers and to waveguides, plain and 
cylindrical, are indicated. New foundations for the design of optical 
instruments, which are centered around an axis, and for the design of 
invisibility cloaks, plain and spherical, are proposed. The variability 
of material properties makes possible effects which cannot happen 
in constant media, e.g. stopping the flux of electromagnetic energy 
without loss. As a consequence, spherical devices can be constructed 
which bind electromagnetic waves. Key words: Electromagnetism; 
Optics; Wave Optics; Diffraction and Scattering; Polarization; Optical 
Materials; Optical Elements and Devices; Fiber Optics 
PACS numbers: 41.; 42.; 42.25.-p; 42.25.Fx; 42.25.Gy; 42.25.Ja; 42.70.- 
a; 42.79.-e; 42.81.-i 

(The final version appeared in Z. Naturforsch. 67a, 111-131 (2012).) 



1 Two steps towards reality 

Let us trace the way Schrodinger walked to find the Schrodinger equation. 
He thought that the propagation of light, if it is construed as propagation of 
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particles, is best described by the eikonal equation: 



(Vs(r))2 = n2(r) . (1) 

Surfaces of equal eikonal s(r) are perpendicular to the light rays everywhere 
in the space described by the vector of location r. n{r) is the index of 
refraction. Schrodinger compared this with the Helmholtz equation 

V^^{r) + !^u;2^(r) = (2) 



c 



which he considered as the best description of waves, u denotes the frequency 
of that wave and c is the velocity of light. The meaning of ip is not known. 
Next Schrodinger remembered that there is an eikonal equation for massive 
particles, too, the Hamilton- Jacobi equation 

(V5(r))2 = 2m{E - V{r)) . (3) 

The mechanical eikonal S{r) has a similar meaning as in ray optics, but its 
dimension is different. Hence Schrodinger deduced from a comparison of ([T]) 
and (Bl) a mechanical index of refraction 



n{r) = ^^2m{E-V{r)) . (4) 

The factor in front of 2m{E — V{r)) is an adjustable constant to get dimen- 
sions right. That it is related to Planck's constant h, Schrodinger realized 
when he solved the first problems. However, when the guess (jl]) is used in 
the Helmholtz equation 02]), an equation for the wavy propagation of massive 
particles is established, the Schrodinger equation: 



Here also the meaning of \E' is not clear. 

The problem with this type of approach is the Helmholtz equation ([2]). 
The wavy propagation of light is reigned by Maxwell's equations. In the 
analytic solution of these equations, the Helmholtz equation occurs as a 
mathematical auxiliary but this is only true when material properties 
as permittivity e and permeability n are constant. What is the replacement 
of the Helmholtz equation if these properties and thus the index of refraction 
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vary in space? This question will be answered in Sections [3] and |H The 
modifications will turn out so severe that the Helmholtz equation can be 
considered only in rare cases as an approximation. It is not even possible to 
formulate the true equations using the index of refraction only. Permittivity 
and permeability enter individually. 

The true companion of Maxwell's equations is Dirac's equation. Maxwell's 
is for vectors, Dirac's for spinors. Yet both systems carry similar informa- 
tion, namely equations for divergences and curls related to time derivatives. 
Dirac's equation is a linear system of partial differential equations with vari- 
able coefficients, the electro dynamic potentials. In Maxwell's equations, vari- 
able coefficients appear when permittivity, permeability and conductivity 
depend on location. Dirac's equation can be solved analytically if the coeffi- 
cients vary just one-dimensionally or if they vary central-symmetrically. The 
analog for Maxwell's equations, and more, is the main result of this article, 
see Sections l3.ll and l3^ 

The impact of this article should be even larger on practical problems. In 
modern times, people fabricate graded materials or so-called metamaterials 
within which permittivity and permeability vary in space almost arbitrarily. 
Therefore analytic solutions that predict effects of such variations will be 
useful. 

Yet usefulness for practicians coerces the consideration of dissipation and 
dispersion. Most materials have finite conductivity. Ohmic currents must 
be included in the theory. Moreover inertia and friction within the materi- 
als modify permittivity e, permeability /z and conductivity a. The simple 
constants must be upgraded, in a minimum approach to reality, to response 
functions which vary in space and describe retardation: 



D(r, t), E(r, t), B(r, t), H(r, t) and j(r, t) denote dielectric displacement, 
electric force field, magnetic force field, magnetic field strength and electric 
current density, respectively. They all are vector fields depending on space r 
and time t. 

With the constitutive relations ([319]) the evolution of the electrodynamic 
field is completely conceived by Maxwell's equations 






(9) 




(10) 
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VB(r,t) = 0, 
VxH(r,t) = 9tD(r,t)+j(r,t) , 
VD(r,t) = p(r,t) , 



(11) 
(12) 

(13) 



written in an unfamiliar sequence for reasons that will become clear in Section 

El 

So these are the two steps to reality: First, solutions of Maxwell's equa- 
tions shall be found with material properties that vary in space. Second, 
retardation shall be taken into account. 

Yet generality will be restricted in two ways: First, only the homogeneous 
problem will be tackled. For example, externally driven currents will be 
omitted. The reason to keep nevertheless the charge density p(r, t) and the 
current density j(r,t) in Maxwell's equations is to admit Ohmic currents. 
The exclusion of nonhomogeneities is not a serious limitation as there are 
standard procedures to construct the solutions of nonhomogeneous equations 
from the solutions of the homogeneous system. 

By contrast, the second lack cannot be cured and can be justified only 
by the desire to produce exact solutions of Maxwell's equations: All material 
properties will be restricted to depend on one or two spatial coordinates 
only. To be specific, introduce coordinates ^, t], ( to describe the vector of 
position r. They may be the Cartesian coordinates x, y, z, but generally these 
greek letters are meant to describe curvilinear yet orthogonal coordinates, for 
example spherical or cylindrical ones. The normalized basis vectors shall be 
denoted as e^, e^, and the line element ds be given as 

{dsf = g^^idO' + 9r,,{dvf + gcddCf (14) 

with elements g^^, gr,rj, fi'cc metric tensor p!8]. The response functions 

are supposed to depend on ( only 

e{C,t), fi{C,t), a(C,t), (15) 

see Section [31 or only on r] and ( 

£{vX,t), f^{vX,t), a{ri,C,t), (16) 

see Section m In the second case (fT6|) . which appears to be more general, 
we will have to impose restrictions upon the solutions of Maxwell's equa- 
tions. Nevertheless the realm of exactly solvable problems will be extended 
immensely. 

This is the plan of this article: In Section [2l Maxwell's and the constitutive 
equations will be rewritten to facilitate a simple description of retardation. 
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In Sections [3] and H] the main results will be produced and proven, namely 
two theorems of representation. They reduce the eight mingled Maxwellian 
equations to two uncoupled partial differential equations each for one un- 
known only. Applications of these theorems are sketched in the Sections 13. 
[3Xn 13X21 KTM K2\ 13X11 lO and US Finally in Section E attempts are 
made to do justice to precursors of the ideas and the results presented here. 

2 Reshaping Maxwell's equations 

The first two equations (11 QUI II) are ideally simple. The aim is to rewrite the 
last two equations f ll21IT3]) until they take the same shape as the first two. 
The clue is to introduce the complete displacement 

C(r,t) = D(r,t)+ rj(r,r)dr . (17) 
Jo 

This takes to 

VxH(r,t) = 9iC(r,t) (18) 
having up to a sign the same structure as ( ITU]) . Similarly (|T3]) appears as 

VC(r,t) = (19) 

having the exactly same structure as (II ip . Because of the continuitity equa- 
tion 

p(r, t) = -V r j(r, r)dr + p(r, 0) (20) 
Jo 

equation (|T9|l is true if there is no initial bunching of charges p(r, 0) = 0. The 
effects of an initial bunching of charges can be covered by a scalar potential 
in a manner explained in [TJ Sec. 2]. It is not extraordinary enough to be 
treated here. 

The reshaped Maxwell equations are now (llOtllip and (I18lll9p . To close 
the system, we have to combine from dTj) and (Q the constitutive equation 
for the complete displacement. It is 

C(r,t)= re(r,t-r)E(r,r)dr (21) 
Jo 

e(r,t) = £(r,t)+ /V(r,r)dr (22) 
Jo 

with the complete permittivity e(r, t). 
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The reshaped Maxwell equations f lTUlfTT]) and (I18til9p together with the 
constitutive equations ([HD and fl?Il) form a closed system, but it is a system 
of integro-differential equations. Performing Laplace transforms 

Ur)= / f{r,t)e'^'dt (23) 
Jo 

where /(r, t) may denote any component of the vector fields or any response 
function, one gets rid of the integrals. The convolution theorem converts the 
integrals in the constitutive relations ([8]) and (12T|) to products: 

B^(r) =/i^(r)H^(r) , (24) 
a(r) =e.(r)E^(r) . (25) 

The Laplace transform of the complete permittivity follows from fl22]) 

euj{r) = e,^(r) + iauj{r)/u . (26) 

It is just a linear combination of the Laplace transforms of the ordinary 
permittivity and the conductivity. 

Also one gets rid of the derivatives with respect to time t: 

POO 

/ dj{r, t)e'^'dt = -luUv) - /(r, 0) . (27) 

The second term on the right-hand side is a valuable peculiarity of the 
Laplace transform as it facilitates straightforward solutions of initial-value 
problems. 

The Maxwell equations fllOHllj) and fll8H19l) are transformed to 



V X E^(r) =zwB,(r) , (28) 

VB,(r) = , (29) 

V X H^(r) = -za;C^(r) , (30) 

Va(r)=0. (31) 



Here two other nonhomogeneities were omitted, viz. B(r, 0) and — C(r, 0) on 
the right-hand sides of (128|) and (l30i) which arise from the Laplace transforms 
of — 9iB(r, t) and 9tC(r,t), respectively, according to equation (1271) . 

Most people would consider the reshaped Maxwell equations fl28H31l) as 
obtained from the orginal Maxwell equations just by separation of exp(— zwt). 
Yet this point of view hides the origin of the permittivity ^^{y) and perme- 
ability /ia;(r) depending on frequency and it aggravates the solution of initial- 
value problems, i.e. it impedes a rational theory of pulses. To do this, one 
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has, first, to solve the reshaped Maxwell equations (128113 ip with the reshaped 
constitutive relations ( 121112^ . second, to introduce the initial- values as non- 
homogeneities in (128|) and (130|) and to solve the nonhomogeneous system and, 
third, to calculate the pulses from the inverse Laplace transform: 

1 roo+ir 

f{r,t) = — Ur)e-^^'du. (32) 

^Vr J —oo+ir 

r denotes a real number big enough such that all locations of sigularities of 
/^(r) in the complex plane of u have smaller real parts. 

The customary variable of Laplace transforms is p = iu. The author 
introduced u instead in order to pacify conservative readers. It they want to 
believe that the equations (|2HHSI]) are just the ordinary Maxwell equations 
with exp{—iujt) separated off, they can do so. The calculations to be pre- 
sented right now, however, do not depend on this point of view. So let us 
abbreviate: 

C = a(r), E = E^(r), e = e^(r), 

B = B^(r), H = H^(r), /x = /.^(r). (33) 



3 Triple curl again 

The aim is to reduce all Maxwell equations f l28ll3T]) to one partial differential 
equation for one scalar auxiliary, the representative b = 6t^(r). The approach 
is the same as in [3 Sec. 2] looking for certain equations with triple curl. 
The ansatz 

B = -V X V X v6 (34) 
E = -iuV X vb (35) 

solves two Maxwell equations immediately, viz. (^E^ and fl2^ . The vector 
field V, the carrier, shall be chosen such that the remaining two equations 
(130|) and (l3Ti) . too, can be solved. Such a choice will be possible if the response 
functions depend only on one spatial variable, say (, as declared in ( IT5|) . and 
thus 

e = ^M), = /^a;(C)- (36) 

Inserting the constitutive relations (I24ti25p into the ansatz (I34ti35p gives 

H = --V X V X v6 , (37) 
C = -iuV X ve6 . (38) 
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In the last equation the permittivity e was drawn under the curl though it 
is not constant. This can be justified if the carrier v is chosen to point into 
the direction of the basis vector e^. 

V = |v|e^-. (39) 

Then, because of fl5B]) . the carrier points into the same direction as the gra- 
dient of e 

„ er de v de , , , 

y/m dC dC 



Consequently in the identity 

eV X v6 = V X ve6 + v6 X Ve (41) 

the last term is zero and thus equation (155]) proven. 

Because of fl38l) . Maxwell's equation f l3T]) is automatically fulfilled, too. 
So we just have to care for fl30l) . Insertion of fl371) and fl38|) produces an 
equation of triple curl 

V X (--V X V X v6 + vew^fo) = . (42) 

The second term behind the leading curl is proportional to the carrier. All 
that remains to be done in order to obtain the desired scalar equation is to 
show that the first term is a gradient plus a term which also aligns with the 
carrier. First, we replace the double curl with the Laplacian V^: 

--Vx Vx v6 = -(VV6-V(Vv6)) . (43) 
/i /i 

Second, we commute the carrier v with the Laplacian and require that the 
commutation does not produce more than a gradient. This can be done only 
if the carrier varies at most linearly 

V = vo + t;ir ^ V\b = vV^6 + 2viVb , (44) 

vq denoting a constant vector and Vi a constant number. The equivalence 
is valid only if the dependence of b on the coordinates is not restricted. For 
details of the proof see [10] or [6]. If fi were constant, we had completed the 
task. Then the second terms on the right-hand sides of (H3|) and (jH]) are 
gradients which the leading curl in fH2l) discards. When fi varies, we must 
effect a third transformation 

-(-V(Vv6) + 2viVb) = (v-) (Vv6 - 2vib) - V-(Vv6 - 2vib) . (45) 
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Here, at last, the second term on the right-hand side is a gradient, while the 
first aligns with the carrier because of fl5B|) . Using again fl5^ we find 

Vi = ^ Ai (46) 



similar to Collecting ([43]), (gl]) and (|45]) we derive from ([42]) the scalar 
equation 

+ --^ — (47-) (Vv6 - 2vib) + tfiio^b = . (47) 



|v|v^VdC/i 

To find formulae for the other polarization, we stipulate the ansatz 

C = V X V X va (48) 
H = -iuV X va (49) 

with the representative a = aaj(r), which automatically satisfies all Maxwell 
equations except (128|1 . The demand to have also this one solved produces a 
triple-curl equation similar to ( H2|) . a replaces b, e and fi are interchanged. 
Performing the same transformations as before, we arrive at a scalar equation 
similar to (H7|) . 

Thus we finished the proof of the 

Three-Dimensional Representation Theorem. Solutions of Maxwell's 
equations / [28lf3l]) are provided by the representations 

E = iv X V X va - iuV x v6 (50) 
e 

H = —iooV X va V x V x v6 (51) 

if the representatives a and b obey the scalar partial differential equations 
V^a - A7(vVloge)(Vva - 2via) + efiu^a = (52) 



- — (v V log /i) ( V v6 - 2t;i6) + eficu^b = (53) 
v^ 



aiid the carrier is chosen such that 

V = vq + fir and Ve oc v and V/i oc v (54) 

vq being a constant vector and vi a constant number. 

In most optical instruments, different media meet. Often a graded medium 
is surrounded by air. Usually material properties jump discontinously where 
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the media touch. Therefore boundary conditions for the representatives are 
necessary. 

Corollary on boundary-value conditions. Let S denote the surface 
where different media meet, n the normal on this surface and d/dn the 
differentiation along this normal. If v oc n, the representatives a and b must 
satisfy 



a\s- — a \s+ , 



ld\\-\a ld\v 



a 



s- e 



dn 



e dn 

ld\v\b l^\^r\b 

b\s- = b\s+, — ^ — = — ^ — 
jj, on yU on 



(55) 

s+ 



(56) 

s+ 



The symbols S— and S+ indicate that the values of the functions and their 
derivatives are to be calculated via an approach on the one side of S, say, 
the low side S—, or on the other side, say, the high side S+. 

Proof: It follows from the structure of the Maxwell equations fl28|) and 
fl5U]) that the tangential components of the magnetic field strength H and 
the electric force field E don't jump on transition through S. Evaluating 
these facts in the representation formulae (15011511) produces the proof of the 
corollary. Details of the calculational procedure are similar as in [Tj Sec. 3]. 
Q.E.D. 

In optics it is difficult to observe the electromagnetic field directly. Instead 
one measures the fiux of energy which can be calculated as the Pointing vector 
S(r, t) = E(r, t) X H(r, t). This is a general formula for physical fields E(r, t) 
and H(r,t). Before we can use it for the mathematical fields handled here, 
we must calculate the dependence on time from fl32|) and extract the real 
parts 3ft. However, if the dependence on time can be described by the factor 
exp{—iut) with real frequency u, we may apply the 

Corollary on the energy flux. The time-averaged Pointing vector S can 
be calculated from the representatives a and b according to 

S = i3ft(E X H*) = 3f^(^(V X V X va) X (V X va* 

+^( —(V X v6) X (V X V X vb* 
V2/i 

+3ftf y(V X v6) X (V X va*) 

X V X va) X (V X V X v6*)) (57) 

the asterisk * denoting complex conjugation. 

The proof follows immediately from the representation formulae (15011511) . 
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When there is only one polarization, i.e. either a = or 6 = 0, the mixed 
terms in the third and the forth lines of flFTI) do not apply. In the first and 
second lines, reader's attention shouldn't miss the inconspicuous imaginary 
units i and the factors 1/e as well as The former are indispensable for a 
weird stop of energy flux in dielectric materials, whereas the latter may alter 
the type of the flux considerably when they are non-constant. 

Due to the condition two limiting cases stand out, namely when 

the material properties vary one-dimensionally, see Section 13. 1[ or central- 
symmetrically, see Section [3^ 



3.1 One-dimensional variations of material properties 

z be the name of the coordinate along which permittivity, permeability and 
conductivity are allowed to vary. It is the same z which is costumary in the 
cartesian system and all cylindrical coordinate systems. In the equations of 
Sections of [T]and|3]we got to set C = The unit vector along z is chosen as 
carrier, i.e. v = e^. According to condition fl54p we have vq = and vi = 0. 
The differential equations ( 1521155]) become 

„9 d log e (9a 9 , , 

V^a -^— + erLU^a = 0, 58 

dz oz 

V%-^^ + e^^u% = 0. (59) 
dz oz 

These differential equations are valid in all coordinate systems which incor- 
porate a cartesian direction. 

The most elementary example is the cartesian system x, y, z. fl58p appears 

as 

d'^a d'^a d'^a dloge^a ^ 

Tr^ + Tr^ + TT^ A — ^ + e/i^a = 0, (60) 

ox^ oy^ oz^ dz oz 

which is separated by the ansatz 

a = XYZa with X = X{x),Y = Y{y), Za = Za^{z) (61) 
to yield three ordinary differential equations 



d^X 



dx"^ 
dV 



+ fc^X = (62) 
2 + fcjr = (63) 



d^Zg dlogedZg ^ ^^^^ ^ ,.2 
dz'^ dz dz 



dy 

+ {eiJUJ^ -e)Za = (64) 
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with separation constants k^, and ky meaning physically wave numbers, and 
with k"^ = k"^ + ky. The solutions of the first two equations fl62H63l) are, of 
course, exponentials or sines and cosines 

a = e'^^'^^'^'^yy^Za , (65) 

but the solutions of (164|) can be weird and convey physical information never 
considered before. 

According to (l59l) . the representative h is subject to another partial dif- 
ferential equation 

d% d'^b d% d\og iidh 2^ g /gg\ 

dx'^ dy'^ dz"^ dz dz 

When this is separated using the ansatz b = XYZ^ similar to fl6ip , the 
ordinary differential equations for X and Y the same as fl62l) and fl63p . re- 
spectively, but the differential equation for Zb 

d^Zfe dlog/i dZj, 2 ,2w n fR^\ 

^-^^ + (^/^- -k)Zb = (67) 

is different from (1M|) . This reflects a responsiveness of graded materials to 
polarization. Examples will be discussed in the next Sections. 

It should never be forgotten that the partial differential equations f l52ll53|) 
hold in any cylindrical coordinate system. For a less trivial example let us 
select elliptic-cylinder coordinates ^, //, C [El Sec.l] 

a; = c cosh ^ cos ?7 (68) 
y = c sinh^ sini] (69) 
z = C (70) 

c being a positive constant. The ordinary differential equations obtained by 
separation of 

a or 6 = SHZa or 6 with S = E{^), H = R^r]), Za or b = Za or 

(71) 

are 

d^- 

-{q-k^c^cosh^O^ = , (72) 



2 



^ + {q- kh^ cos^ r/)H = . (73) 
dr]^ 
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The separation constants are here q and k^. The solutions of these two 
equations are Mathieu functions [T7]. The third equation was not written 
because it is identical with flMj) or ( 167|) . The elliptic-cylinder coordinates 
are especially interesting as they allow to exactly predict the diffraction by 
a strip or a slit, not just by an edge, see [3 Sees. 10, 11] and the references 
therein. Moreover in the theory presented here, the strip need not to be 
homogeneous. So we can devise novel ways to bunching and debunching of 
electromagnetic waves. 



3.1.1 Graded reflectors, transmitters and polarizers 

As a first example of application, consider a medium homogenous and isotropic 
for 2; < 0: 

e = e_ , = /i_ (74) 

with constant permittivity e_ and permeability //_. The solutions of the 
differential equations (160|) and (166|) are almost trivial: 



a_ or 5_ = e'^>^Myy)^^i^7ii^::fl^z + ^ ^-i^/Ti;:::^ ^ ^^5^ 

They describe waves of different polarization coming from negative infinity 
and being partially reflected at z = 0. Ra and Rh are complex constants to 
fix the strengths and phases of the reflected waves. These constants have to 
be determined from the solution of the boundary-value problem; see below. 
For z > Q the permittivity is supposed to vary 

^i = fi+ , (76) 



1 - (1 - eoo/e+)exp(-2;/z+) 

whereas the permeability is supposed to stay constant The permittivity 
fl76|) takes the value e+ at z = and increases or decreases with positive slope 
constant z+ while it approaches eoo for 2; — t- oo. 

The ordinary differential equations flM]) and fl57|) differ considerably. Yet 
both can be solved in terms of the hypergeometric function 

or 6 — e V 

F[pLa or fe,/3a or 6,7a or b, (1 - eoo/e+)exp(-2;/z+)) . (77) 
For Za we have to use the parameters 



= 1/2 (1 + a/1 + 4^24) - Veoo/i+c^2 _ A;2 



Pa = 1/2 (1 - ^1 + 4P4) - i^t^iJi+uj^ -k'z+, (7J 



7a = 1 - 2z JeooA*+w2 - k'^ z+ , 
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but for 



(79) 



The very fact that the parameters in flTHl) and (1791) differ shows that graded 
materials act discriminatorily towards waves with different polarization. 

There are second solutions of the second-degree equations ( 1M|) and (IB7|) . 
but we don't need them here. Namely 



F(a,/3,7, (1 - eoo/e+) exp(-2;/2;+)) 1 for z 00 



^0) 



is a property of the hypergeometric function for all values of a, /3, 7. Therefore 
(177|) is sufficient to describe a wave running to positive infinity. 

Thus we have, also for z > 0, solutions of the differential equations ( BUj) 
and (1661): 



a_ or 6_ = e*(*^-^+'=^^) T„ or 

i^(aa or fe,/3a or fe,7a or fe, (1 - eoo/e+) exp(-2;/2;+)) . (81) 

Ta and Tb are complex constants to seize the strengths and phases of the 
transmitted wave. 

The four coefficients Ra,Ta and Rb,Tb are determined by the conditions 
at the boundary S. The surface 5* is given here as z = 0. The normal n 
coincides with the vector e^. Specializing f lS^llSB]) yields 



Za \z=-0 '- 
Zb \z=-0 



Za \z=+0 ; 
■ Zb \z=+0 ; 



IdZa 


_ IdZa 




e dz 


z=~o e dz 


z=+0 


IdZb 


_ IdZb 




fi dz 


z=-o jJ' dz 


z=+0 



^2) 
^3) 



The boundary conditions fl82p and flS^ yield both two linear equations 
for Ra,Ta and Rb,Tb, respectively, and are solvable elementarily. Next one 
differentiates the electromagnetic field from the representatives ( 175|) and ( IHTj) 
according to ( l50ll5T|) and calculates from the electromagnetic field the Point- 
ing vector to obtain quantities directly comparable to experimental results. 

Several modifications and generalizations are at hand. The waves may 
run in the opposite directions. In this case the second solutions of (1M|) and 
fl67|l are needed to describe reflected waves. They too are expressible in 
terms of the hypergeometric function. We may also study the transition of 
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waves from one graded medium to other graded materia. In this case the 
purely exponential waves f l75p must be replaced with expressions containing 
hypergeometric functions. 

In the definition of gradation (1761) . it was assumed that only the per- 
mittivity varies. This is resonable for many materials, but for independent 
variations both of permittivity and permeability it is advisable to use other 
parametrizations. Power laws, for example, 

e = e„;z", ^ = ^ipz^ with a + /3 = -2, -1, 0, 2 (84) 

lead to solvable differential equations (IMj) and (^7\i - solvable by functions 
not more complicated than the confluent hypergeometric function. 
Another useful variability is 

z z 
e = Co exp — , fj, = fiQ exp — (85) 

Zf Zf, 



with positive constants eo, /io and positive or negative constant decay lengths 
z^ and z^. Apt solutions of (|M|) and 
of weird index and weird argument: 



z^ and z^. Apt solutions of (IM|) and (16 7p are Hankel functions H^^ and H^^ 



z \ ^(1.2) f '2zeZ^,^/e^u: {z, + z^)z 



z^Jl + Ak^z^ 



with u = , (86) 



Ze + Z^ 



Z = exp f —] ^2^.^/.v'e^a; _ {z^ + z^)z 



b = exp ( — I ii;; ' M J exp 

^Z^ / \ -r Z^ ZiZf^Z^ 



with V = . (87) 

Ze + Zf, 

Hankel functions, also denoted as Bessel functions of third kind fTEl Chap. Ill], 
are more appropriate for the description of running waves than ordinary 
Bessel and Neumann functions. 

The parametrizations (18^ or (185|) should be applied to slabs as there is no 
medium with infinitely large or infinitely small material properties. Instead 
of one block of boundary- value conditions as in (I82ti83p there should be two 
blocks, one for a boundary at, say, z = Zi and the other block for a surface 
at ^ = ^2- The somewhat larger linear systems do not essentially aggravate 
the solution of the entire problem. 
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3.1.2 Stopping the energy flux 

A peculiar special case of (185|) is 

z — z 

e = eo exp — , = /iq exp — . (88) 

Ze Ze 

The index of refraction ([6]) is just 1. Believers in the Helmholtz equation ([2]) 
should expect just ordinary propagation of waves, but no spectacular effect. 
The true differential equations flM|) and f E7|) have constant coefficients: 

^ - + (eo/xoc^^ - k')Za = , (89) 

dz^ z^ dz 

Hence their solutions are readily found: 



Za = exp (^^^ exp ± y ^ - (eo/iow2 - k^) z^ , (91) 

Zb = exp (^^^ exp ± ]Jj^- (eoyUoc^^ - /^^) 2;^ . (92) 

The exponential functions directly behind the equal signs impress much, but 
for the energy flux they don't matter at all. They only exist to compensate 
the factors 1/e and l//i in ( 157|) . The astonishing item is the square root 
in the second exponentials. Normally, i.e. for z^ ^ 00 as it holds for any 
almost homogeneous medium, the value of the root is imaginary and the 
exponential function carrying it is complex, just the description of a plane 
wave propagating we are used to. Yet if 

l^^l < o / \ ^2 ' (93) 

the functions (19T]) and (192|) become real and the energy flux against the 
gradation is stopped. The condition means that the slope constant must 
be smaller than the wavelength the wave would have in a medium without 
gradation divided by An. Though this is a short length, it can be constructed 
in modern labs. 

Proof: Evaluation of fl57|) using fl9T|) and fl65l) yields 



S oc e^A;^. + e^A;^ + e, <! ^Jeofiou^ - - 1/ {2z,)^ if eo/io^^ > fc^ + l/(22;,)2 

otherwise 

(94) 
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The result is the same when the electromagnetic wave is represented by b 
with Zb from (1^. 

The most surprising feature of the effect is its independence of the sign of 
the slope constant z^. It doesn't matter if permittivity increases or decreases. 
Important is only a sufficiently steep change. 

The finding fl94p differs fundamentally from the behavior of electromag- 
netic waves in conducting materials. There the energy intrudes and is dis- 
sipated. The finding is also fundamentally different from the behavior in 
homogeneous dielectric materials. Even when there is a discontinuity, part 
of the wave is maybe reflected, but the remainder goes on to transport energy. 
The waves found here differ as well from the evanescent waves which make 
possible dielectric waveguides. In these waveguides, k must be greater than 
a positive cut-off wavenumber, whereas the stopping described here works 
also for k = 0. 

The effect is moreover not singular. The case (l88i) need not be fulfilled 
exactly. One can derive this from the asymptotic expansions of the Hankel 
functions in ( 15B1I57I) . When both index and argument get great, the Han- 
kel functions become exponential functions with real argument if the index 
is greater than the argument, but they become exponential functions with 
imaginary argument in the opposite case [T6l Sec. 3. 14.2], similar to the ele- 
mentary functions in fl9HI92p . Generally, however, there is some dependence 
of the stopping on polarization. 

Time-dependent analysis reveals that the energy density oscillates. Dur- 
ing one half of the period 27i/u it is pushed into the graded medium, during 
the other half it is withdrawn. The depth of the penetration is approximately 
described by the second exponential functions in (19111921) taken with negative 
signs before the roots. Therefore, in a slab of finite thickness, the stopping is 
not perfect. Waves impinging on the one boundary of the graded medium de- 
crease in the medium, but the mechanism just described may excite waves, 
though weak ones, on the other boundary. Exact amplitudes and phases 
follow from the boundary conditions fl82ti83p . 

3.1.3 Dielectric mirrors 

Having read the previous Section one may argue that monotonous exponen- 
tial growth cannot be maintained on long distances. However, one can realize 
similar stopping with zigzagging material properties. This section is devoted 
to periodic variations of the permittivity. Let 

e = eo — ei cos koz , /x = /iq (95) 
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with positive constants eo, ei, Hq and /cq- When we compare the resulting 
equation f E7|) 

+ ieQiiQUJ^ - k"^ - ei/iow^ cos koz)Zb = (96) 
to the equation of the Mathieu functions me^{w; q) 



(A — 2q cos 2w)me±y{w] g) = (97) 



with constant A and q as defined by Meixner and Schafke pTf p. 105], see 
p. 404] for a shghtly different definition, we find 

Zb = me±u{koz/2; q), A = 4:{eof^o<^^ - k^)/ko, q = 2eifioUj^/kl . (98) 

The Mathieu functions possess properties transcending the flexibihty of 
the hypergeometric function and all its descendants. It follows from Floquet's 
theorem [251 p. 412] that one can compute them from a Fourier series times 
an exponential factor 

n=+oo 

me,(^;g)=e^^'" ^ c^„(g)e^'"" ■ (99) 

n=—OD 

with coefficients C2n{q) for which Hill's theory [2^ p. 413] provides handy 
expressions. 

The characteristic exponent u is a. surprising function of A and q. Imagine 
g > fixed while A varies. Then u assumes, for certain bands of A, only real 
values. This is what one should expect for physical reasons: A periodic per- 
turbation in the differential equation causes periodic or rather quasi-periodic 
perturbations in the solutions. However, in the complements of these bands, 
u acquires complex values. The phenomenon is known as parametric amplifi- 
cation, but it is often forgotten that only one solution describes amplification, 
whereas the other describes attenuation. This type of damping, which comes 
about without friction or spatial dissipation, is the effect we want to consider 
here. 

Floquet's theorem and Hill's theory were reinvented and generalized in 
solid-state physics where the theory is known as Bloch's theorem|3],[l3]. In 
solid-state physics, A plays the role of the energy in Schrodinger's equation 
and the complementary bands are denoted as forbidden. 

The lowest forbidden band of Mathieu's equation 0971) is characterized 
m P-120] by 

\X - 1\ < \q\ + 0{q^) . (100) 
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Maximum attenuation, described by the imaginary part of the character- 
istic exponent z/, takes place in its median [T71 p. 165] 

3z/= |g|/2 + 0(g^) where A = 1 + O(g^) . (101) 

To estimate the attenuation of the Pointing vector (1571) . we must take the 
absolute square of the leading factor on the right-hand side of fl99|) . The 
attenuation of energy flux is thus 

S oc exp(-|g|w + 0(g2)) . (102) 

Applying this to the problem at hand (198|) . we see the first forbidden 
band approximately defined by 

|4(eo/iow^ -k'^)-kl\< 2eifioUJ^ . (103) 

The median of the forbidden band is where the left-hand side is zero. One 
may write this as 

11 . ^ 

— ^ — - (104) 

ko 2A/eo/io^J^-fc^ 
meaning that the wave length of the intruding light must be twice the wave 
length of the dielectric zigzagging. In other words: We must have four layers 
of different media or two periods of permittivity for every spatial period of 
light. The reader is encouraged to compare (11041) with (jHSD- 

When fll04p is used in the formula for q ( l98l) and if it is assumed, just for 
simplicity, that k = 0, i.e. the light impinges vertically, we obtain as a crude 
estimate 

q^^ (105) 
2eo 



and thus for the attenuation 



Socexp(-^^). (106) 

A variation of permittivity ei ~ 0.13eo appears to be realistic. Then, accord- 
ing to f ll06p . it takes less than 6 periods of permittivity, or 11 layers, to attain 
attenuation by a decade, it takes less than 23 layers to attain attenuation 
by two decades and so forth. An electromagnetic wave that impinges on a 
periodic structure cannot penetrate. It is reflected. 

The dielectric mirror just described is selective. Light that has not the 
suitable wave length fll04p passes. Yet selectiveness isn't overly sharp. Ac- 
cording to fllOOp . the width of the band is 2|g|. If g is estimated according 
to f ll05p . we find that the relative width of the forbidden band is ei/eo, i.e. 
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25% in the present example. This is enough to cover a considerable part of 
the spectrum visible to human eyes. Moreover the dielectric mirror alters 
its properties with the angle of incidence. The angle is contained in the 
transverse wave number k and enters the theory via the parameter A in (l98ll . 

It remains to check the dielectric mirror for its dependence on polariza- 
tion. To this end the solution of the differential equation has to be 
compared with the solution of ( E7^ which we just discussed. 

The term with the first derivatives in f l6^ can be eliminated introducing 
the auxiliary function Za, yet at the cost of more complications in the factor 
of Za. 

Za = Za/ \/e , 

^+1^^" +2^di^-4^ldij j^" = °- ^'''^ 
Generally this is not exactly a Mathieu equation, but it is, because of its 
periodic coefficient, of Hill's type. It can be solved in the same way as Math- 
ieu's and exhibits the same features, namely allowed and forbidded bands. 
Nevertheless for small oscillations of the permittivity ei -C eo, equation fll07p 
can be approximated by the Mathieu equation 

d^Za ( ^\k^\ 

+ (eo/UoW^ -k^ - (^ei/iow^ - ^yj cos koz)Za = , (108) 

i.e. Za is represented by a Mathieu function, too, where A is same as in 
but 

2eiHoUj'^ ei -ei 







(109) 



Repeating the same deliberations as those following fl98|) . we find that the 
median of the forbidden band is at the same position f ll03p . The parame- 
ter q has now, apart from its sign, approximately the same magnitude as 
in fllOSp . but the sign of q can be compensated by an unimportant phase 
shift of the argument in Mathieu's function. Therefore both the width of 
the forbidden band and the attenuation are approximately the same as for 
Zf,. Thus, surprisingly enough, the dielectric mirror depends but weakly on 
polarization. 



3.2 Central- symmetric variations of material proper- 
ties 

For this case we specialize the condition fl5^ by vq = and fi = 1. The 
carrier is just the vector of position v = r and hence = r^. In all equa- 
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tions of Section [2] we got to replace C with the r customary in the spherical 
coordinate system r, 6, ip. The differential equations (15211531) appear as 

Id'^ra ^ f ^ ^ ( ' O^'^^ ^ dloge Idra ^ g fl 

r dr'^ r'^\sai6d6\ dO J sin^ 9 dip'^ / dr r dr 
Id^rb 1/1 df db\ 1 d%\ dlogfildrb 2, 



r dr'^ r'^\sm6d6\ 06 J sin 6 dip"^ J dr r dr 
Both equations can be separated by similar ansatzes 

a = ^Ra Yim and b = ^Rb Yim with 
Ra = RaM. Rb = RbUr), Yi^ = YUe.if). (112) 

The equation for the angular factor Yi^ is the same for both representatives 
a and b: 

-^^(sin^^) + + + 1)^- = • (113) 

sm 9 dO\ 06 J sm^ 9 d(p^ 

It is the differential equation of the familiar spherical harmonics. They must 
be unique. So / and m must be integers, in fact / = 1, 2, 3, and |m| < I. I = 
is excluded because Fqo is a constant, and a representative not depending 
at all on the angles is annihilated by the curls in (150|) or (l5Ti) . 

The differential equations for the radial parts are extraordinary: 

d^Ra dloge di?„/ 2 i([+i)\p n ^^^A^ 

-r^ 1 2 — jRa = 0, 114 

dr^ dr dr V r^ ) 

d^Rb dhgfidRb f 2 i([+i)\p n nip;^ 

-— — + e/iw ^ — )Rb = 0. 115 

dr^ dr dr V r^ J 

Although they look like their one-dimensional analogs in ( 164|) and ( 167|) . the 
last terms on the left-hand sides are different. Nevertheless also the present 
differential equations can be solved for all power laws conforming to 

e = e„r", fi = fi^r^ with a + /3 = -2, -1, 0, 2 (116) 

i.e. solved in terms of functions not more complicated than the confluent 
hypergeometric function. 

3.2.1 Bound electromagnetic waves 

In Schrodinger's quantum mechanics, electrons can be bound in a spheri- 
cal well. The eigenvalues of Schrodinger's solutions are discrete. A similar 
construction for photons isn't known. For instance 

eo if < r < ri, u = / '"^ if < r < ri, 

ei if ri < r < 00. \ /ii if ri < r < 00. 
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with positive constants eo,/io inside the spherical core r < ri and other 
positive constants ei,/ii outside admits as solution of f lll4p only 



-Ra 

r 



Aii{y/e^ur) if < r < ri, 

B h\^\y/eijli Lor) + C hf\yjeiij,i ur) if ri < r < oo. 



:ii8i 



with spherical Bessel functions ji, spherical Hankel functions hl^\h\^\ I = 
1, 2, 3 . . . [21 Sec. 10], and constants A, B and C . All functions describe run- 
ning waves. One might think of materials where the real permittivity or the 
permeability are negative. Such materials exist, but those negative values 
take place only in narrow bands of u and come always with considerable 
conductivity. So there is not the least chance to establish bound electro- 
magnetic waves and discrete values of u with spatially constant material 
properties. 

With graded materials, however, we can construct a home of bound waves. 
Consider instead of f lll7p 



Co if < r < ri, 

eiTi/r if ri < r < oo. 



The respective solution of flll4p is 



Ra 



fiQ if < r < ri, 

yUiTi/r if ri < r < oo. 



if < r < ri. 



:ii9) 



(120) 



i/.ia; r, ^ ^ < r < OO. 

The solution on the flank of the wall ri < r can be a wave, though a weird 
one 

r±^/W^^F^^J^^l = exp{±i^jeiiiiuj^rl - l{l + 1) logr) (121) 
if uj is sufficiently high. Yet if 



y/e^uri < J^^lil + 1) 
V ei/ii 



:i22i 



with VL as nondimensional substitute of w, the flank function in fll20p just de- 
creases without any variation of phase such that the same considerations ap- 
ply as in Section 13.1.21 The transfer of energy through the flanks is stopped. 
One can check this explicitly evaluating the Pointing vector f l57|) with fll20p 
and ffTT2|) . 

To flnd the eigenvalues of VL and thus of w, we must satisfy the boundary 
conditions The surface 5* is now the sphere r = ri and the normal 

vector is n = r/r. Hence 



Ra 



r=r\ —0 



Ra 



r=ri+0 1 



e dr 



r=ri— 



e dr 



r=ri+0 



(123) 
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This produces a homogenous hnear system for A and B. It has a non-trivial 
solution if its determinant is zero: 



co/^i /co/^o 



mm 



(124) 



the prime indicating differentiation with respect to the argument f2. The 
function on the left-hand side is a parabola open to the left. We need its 
positive branch at positive values of Q. The parabola disappears for Q greater 
than the cut-off given on the right-hand side of fll22p . The function on the 
ride-hand side of (11241) takes the value — (/ -|- 1) at = 0. It increases with 
Q and crosses the fi-axis at the zero of The function continues 

to increase until it approaches its pole at the zero of To secure the 

existence of a solution of (I124p . it would be sufficient to demand that the 
cut-off in (I122p be greater than the zero of (Qji^Q))'. But these zeros are 
not tabulated. So let us be generous and demand that the cut-off be greater 
than the first zero of This yields a condition 



which warrants the existence of at least one positive solution of (I124p . The 
expression on the right-hand side tends to 1 as Z tends to infinity. Therefore 
the restriction on the ratio of the indices of refraction is unimportant at high 
multipolarities. Yet even for / = 1 the condition (11250 can be fulfilled. The 
first zero of the first spherical Bessel function is fli ~ 4.5 Sec. 10] such that 
ratios of the indices must be greater than 3.2. In modern times where indices 
of refraction can be made as big as 38.6 [8], this is moderate requirement. 

Waves of the other polarization can found by replacing the representative 
a with b and by an interchange of e and fi. The characteristic equation of 
this case differs from (I124p just by a different leading factor. Therefore it 
depends on the polarization whether an electromagnetic wave can be bound, 
but it does not depend much. 

One might compare the construction explained here with a hydrogen 
atom. Rather it is similar to a nucleon bound in a collective nuclear po- 
tential as the spectrum of eigenvalues is finite. The essential difference, how- 
ever, is that every electron always carries the same charge which necessitates 
a normalization of its wave. Here, by contrast, the energy of the bound 
electromagnetic wave is arbitrary. The only necessity to confine the energy 
is a possible breakdown of material properties ([319]). It is therefore blame- 
worthy to speak about "photonics", an "atom for photons" and so forth. 




where = 



(125) 
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Nevertheless, if quantum electrodynamics were true, the energy stored in the 
construction just described should be discrete. 

As in Section [3. 1.21 people might argue that the system just constructed 
is not realistic. Fortunately there is no singularity at the origin at r = 0, 
but it is certainly questionable to require permittivity and permeability ap- 
proaching zero as in (11191) . The solution of this problem, however, is known. 
One must replace the monotonous decrease with zigzagging as explained in 
Section 13.1.31 This will work. For the differential equations ( 16^ and flll4p 
are the same for r — )■ oo. 



4 An alternative theorem of representation 

In the study of graded fibers, one cannot use the theorem of representa- 
tion provided in Section [31 It is possible to analyze electromagnetic fields 
in cylindrical bodies, but the material properties must not vary except in 
the direction of the axis of the cylinder. For graded waveguides, one needs 
permittivity and permeability varying with the distance from the axis, i.e. 
e = ^uj{p), = f^uj{p) in circular cylindrical coordinates p,f,z. Yet there 
is no carrier according to ( 154]) that would be proportional to Ve and V/x. 
Fortunately we can rely on the 

Two-Dimensional Representation Theorem. In a system of orthogonal 
coordinates C,,ti,( where the elements of the metric tensor 

9^^ = miv, C), 9vv = 9vvi^^ C), 9^ = 9cdv, C), (126) 



cf. the line element i fMj) . do not depend on the distinguished coordinate ^ 
and where permittivity and permeability 

e = ea;(^?,C), P = Pu{r]X) (127) 

do not depend on ^, the fields 

^ 1^ ah , . 

E = -V X e. luoe. (128) 

H = -iuoe.— -V X ej— ^ (129) 



solve Maxwell's equations i\2M31\) including the constitutive relations ( 1241125)) 
if the representatives a and b 

a = a^{riX), h = h^{riX) (130) 
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do not depend on ^ and obey the differential equations 

d_ I g^^ da 



9ii 



d_ 



9riV 



da 



9a9(i; dC 



1 d log e da 



9m 



9ii 



drj drj 

- 1- 

dv V 9m9ii 9r] 



dlogeda ^ 

+ ejj,uj a = U 



db 



9m9a 
1 d log /i dh 



+ 



d_ 



I 9m db 



1 9 log /i db 

g.ny dr] drj g^^ dC dC 



+ efiuj% = . 



(132) 



For the proof, let us start with the representative a only. The ansatz 



C = 
H 



V X eg 



a 



9ii^ 
a 



(133) 
(134) 



is the special case of fll281ll29|) with 6 = 0. It solves at once the Maxwell 
equations f l30ll3T|) . Yet under the geometrical restrictions fll26p . (1127^ and 
fll30p . it also solves Maxwell's equation fl29]) . namely 



VB 







(135) 



because there is nothing behind the differentiation depending on ^. 

Thus the only Maxwell equation that still expects solution is f l28|) . Using 
the ansatz (113311134]) it is transformed to 



eV X -V X eg 



efiu eg 



(136) 



which is equivalent, as we will see soon, to the differential equation (I13ip . 
It is apparent that the right-hand side of fll36p is proportional to eg/y^^. 
We will discover that the same is true for the left-hand side. To this end we 
sever the differentiation of e in (11361) : 



eV X -V X e 



— (Vloge) X V X eg h V x V x eg 



The nabla operator applied to log e yields 



Vloge 



d log e 



91oge 



9a dC 



(137) 



:i38l 
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The simple curl in fll37p is 



V X e,^ = _ _^^da ^^^^^ 



Evaluating the cross product of (11381) and fll39p as required in fll37p produces 
the second contribution to (11311) times / ^/g^- The double curl on the right- 
hand side of (I137P gives the first contribution to (I13ip times e^/y^^. So it is 
shown that solution of (I13ip completes the solution of Maxwell's equations. 
The truth of (I132p can be proven when we start from the ansatz (11281 



I129p with the representative h only, putting a = 0. All Maxwell equations 
turn out to be automatically solved except fl30|l . This one is evaluated as 
described in (I136m39p where a is interchanged with h and e with /x. 

Finally we remember the linearity of Maxwell's equations. The full proof 
of the alternative theorem of representation is just the superposition of the 
two proofs produced in the last paragraphs. Q.E.D. 

The consistent setup of boundary- value problems is described in the 
Corollary on boundary-value conditions. Let S denote the line where 
different media meet, n the normal on this line with ne^ = and d/dn the 
differentiation along this normal. The representatives a and b must satisfy 

1 da 
e dn 



a\s- = a \s+ , 

e dn 



1 da 
e dn 
1 db 

\s+ , --^ 
jjL dn 



1 dh 
/i dn 



s+ 



(140) 
(141) 



The symbols S— and S+ indicate that the values of the functions and their 
derivatives are to be calculated via an approach on the one side of S, say, 
the low side S—, or on the other side, say, the high side S+. 

The proof is nearly the same as the proof of the corollary on boundary- 
value conditions in Section |3l The main difference is: We must use now the 
representation formulae (I128tll29]l . 

The reader is kindly asked not to misunderstand the denotation "two- 
dimensional" . A propagating electromagnetic field always spans the three- 
dimensional space. "Two-dimensional" means just that all components of 
the field depend only on two coordinates. 

Yet quite a few problems can be declared to be two-dimensional by a 
judicious choice of coordinates. In all these cases it is advantageous to apply 
the two-dimensional representation theorem. Namely the calculation of the 
electromagnetic field using (I128p and (I129p takes less work than using ( I5U]) 
and (1^ as two curls less need to be computed. 

When permittivity and permeability are constant in space, the two-dimensional 
theorem of representation doesn't offer anything which is not included in the 
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three-dimensional theorem given in Section [31 With constant material prop- 
erties, the theorem already presented in [71 Sec. 10] grants the best systematic 
approach. 

4.1 Examples in a plane 

The most straightforward applications of the foregoing theorem take place 
in cartesian coordinates x,y, z. None of the components of the metric tensor 
depends on any coordinate: 

gxx = 1, 9yy = 1, 9zz = l ■ (142) 

As distinguished coordinate ^ we may select either x or y ot z. Let us identify 
^ = z,rj = xX = y- The partial differential equations fll3Hll32p appear as 

d^a d'^a 9 log e (9a 
\ 2 

Qy2 Qj. 

d% d'^b dlogfidb 
dx"^ dy"^ dx dx 

The electromagnetic field is represented according to (1128111291) through 

E = -V X e^a — iujej) (145) 
H = —iue^a — —V x e^b (146) 

showing that the electromagnetic wave extends in three dimensions while the 
representatives a = a^{x,y) and b = b^{x,y) depend on two coordinates x 
and y. Moreover it should be noticed that both permittivity and permeability 
may depend on both coordinates x and y. Therefore for the so-called two- 
dimensional problem, the equations f l 143111 441) constitute the most general 
reduction of the coupled Maxwellian system to two uncoupled equations. 

Even when the partial differential equations (1143111441) are not separable, 
they vastly simplify the solution of Maxwell's equations as all methods which 
people learn in the ordinary courses of quantum mechanics can be applied 
directly, for example, Born's approximation and the JWKB, denoted also 
as semi-classical approximation. In the latter case, however, the classical 
eikonal equation will turn out to be only of restricted usefulness. Also 
numerical methods will profit from the reduction. 

When e and fi depend only on one coordinate, x or y, the partial differen- 
tial equations (I143til44p can be separated and produce ordinary differential 



(9 log e 9a 2 n 

— — — + efiu a = , 
oy oy 

d log /i db 



dy dy 



+ eij,u% = . 



(143) 
(144) 
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equations similar to f lM|) and f E7|) . Also when e and fi are products of func- 
tions which depend on one coordinate only, separation is possible, but only 
under certain circumstances. We will see an example below. 

In cartesian coordinates we can choose x, y or z as the distinguished 
coordinate ^, but a changed choice does not alter the geometrical situation. 
In circular cylindrical coordinates p,(p,z the elements of the metric tensor 
depend neither on z nor on (f: 

9pp = 1, 9w = P^ 9zz = 1- (147) 

Hence we may select either or 2; as distinguished coordinate, but now the 
choice varies the geometrical situation. 

Begin with ^ = z as distinguished coordinate. The partial differential 
equation f ll3ip appears as 

1/9 da\ 1 d'^a dlogeda 1 dlogeda 2 g (14,S) 

pxdp^dp) p'^ d(f'^ dp dp dip dip ^f^^ ^ 

The equation fll32p for b is up to an interchange of e with p identical and is 
therefore not written. 

The partial differential equation fll48p is the equation of plane scattering. 
It can be separated if e and p are functions of p only. The ansatz 

a = P„ $ with P„ = FaM , ^ = Hv) (149) 
generates the ordinary differential equations 



d^$ 



dip 

d^Pa fl dloge\dPa / 2 



+ m2$ = (150) 



+ -^+{epu'-— Fa = (151 



dp^ Vp dp / dp V p 

The first equation has the familiar solutions $ = exp{im(p). The separa- 
tion constant must be the square of an integer m = 0, ±1, ±2, ±3, . . .. 
Otherwise a and hence the electromagnetic field would not be unique. The 
second equation can be solved using functions not more complicated than 
the confluent hypergeometric function if 

e = e«p", p = ppp^ with a + /3 = -2, -1, 0, 2 (152) 

and with constant a, /9, and pp. This together with the boundary condi- 
tions (11401) gives ample freedom to model graded centers of scattering around 
p = 0. 
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Interestingly the partial differential equation f ll48p can also be separated 
if permittivity and permeability depend on both variables, e.g. 

e = p"e^, fi = p^fi^ with a + = -2 and = e^^{(p), p^ = p^u{<^). (153) 

The functions e and p must be periodic for the uniqueness, but else they are 
arbitrary. The ansatz 



produces 



a = Fa^a with P, = P„(p), = <^Uv) (154) 

d^^a dloge^d$ 2 
1^ - + + W = , (155) 

d^Pa 1 - a dP„ g , , 

dp^ p dp p"' 

with the separation constant q. The latter equation is now of Eulerian type 
simply solved by 

P, = p"/2±V'?+°V4 ^ (157) 

whereas the former equation is of Hill's type. When the average values of 
e,„ and denoted by and p<^, respectively, and if and p^ oscillate 

around their average values but weakly, then fll55p is solved by periodic 
Mathieu functions and the separation constant can be estimated as 

- e^p^Lo^ (158) 

m = 0, ±1, ±2, ±3, . . . Because of (11571) this determines whether we see a 
wave or monotonous variation along p. 

There is no sizeable difficulty to solve f ll56p even if the oscillations of 
and large. Use, for example. Hill's theory. 

Presently the devices to simulate invisibility cloaks are mostly plain. 
Therefore this is the Section with the best formulas to design them. The 
best formulas to design central-symmetrical cloaks can be found in Section 

4.2 Examples around an axis 

When we stay with circular cylindrical coordinates, but select ^ = cp as 
distinguished coordinate, we obtain from (1131111321) 

/d lda\ ^ d^a dloge da dloge da ^ ^ ^2^ = (159) 

^Kdppdp) dz"^ dp dp dz dz ^^'^ 

/d ldb\ ^ d% dlogp db dlogp db ^ ^ ^% = q ngg) 

^Xdppdp) &2 dp dp dz dz 
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The reader might notice that the second order operator in these equations 
cannot be understood as a part of the Laplace operator. Also one needs to 
get used to the equations of representation 

E=lf-^|^ + ^|^V^^-& (161) 
e V p oz pOpJ p 

1 ( e.p dh dh \ 



H = -tuj^ a+-(^^-^^] (162) 
p p \ p Oz p Op J 

which follow, despite of their weird appearance, directly from f ll28til29]) . For 
6 = the magnetic field forms rings around the axis 2 = 0, for a = it is 
the electric field which clings to circular lines. 

Most optical instruments are centered around an axis. The equations 
fll59m62]) provide a better foundation to design them than anything known 
up to now. First, their solutions yield exact solutions of Maxwell's equations. 
Second, one has to solve partial differential equations only for one unknown. 
This is much simpler than solving the multiply coupled Maxwell equations. 
One may, for example, insert a singularity a or 6 ~ log((p — po)^ + (2; — -Zo)^ at 
a point of an object po, Zo and compute where other singularities arise. One 
may simulate lenses by regions in the p — z plane where permittivity e and 
permeability p are increased. At the same time, one may similate metallic 
stops of finite thickness by regions with complex e, cf. ( 126|) . and study the 
interaction between lenses and stops. Many properties of imaging can thus 
be predicted, however, with an important exception: astigmatism can not be 
observed because the dependence on the azimuth is missing. 

The first example is the graded mono-mode fiber. There we have permit- 
tivity e = etj(p) and permeability p = Puj{p) as functions of the axial distance 
p only. Both partial differential equations f ll59p and fll60p can be separated. 
Let us select the first for example. The ansatz with leading p 

a = pFaZ with = P,,(p), Z = Z{z) (163) 

is advantageous because it permits simple boundary conditions for p — 
and p — 7- 00, namely Pq = in both cases. Why? We just have to consider 
the electromagnetic field in fll6mi62l) and to demand that all its components 
should stay finite on the axis and decrease towards infinity. The ansatz fll63p 
generates these ordinary differential equations: 

+ k^Z = (164) 



d^2 

+ {--^^]^+{epoj'-k'-^--^^]P, = (165) 



d'P, , fl dloge\dPa / 2 ,2 1 Idloge 



dp^ Vp dp / dp V ^ p^ p dp 
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The relation between the real wavenumber k and the frequency 00, the so- 
called dispersion relation, is the desired item. It is found from the solution 
of f ll65p satisfying the boundary conditions for p — )■ and for p — > 00 as 
explained in the previous paragraph. One can solve the ordinary differential 
equation fll65p for permittivities and permeabilities obeying power laws as 
in (11521) or one can solve it numerically. A boundary-value problem with 
one ordinary differential equation is by orders of magnitude simpler than the 
same problem with the Maxwell equations and it is much simpler to attain 
high accuracy. 

The considerations produced here for circular cylindrical coordinates can 
be transferred to all coordinate system which embody an axis of rotation, e.g. 
spherical coordinates, prolate and oblate spheroidal coordinates, parabolic 
coordinates and all rotational systems [ISl sees. I, IV]. Especially interesting 
appear at first glance the oblate spheriodal coordinates because they allow 
an easy study what an electromagnetic wave does in a bottleneck, and the 
toroidal coordinates because they allow the study of electromagnetic waves 
in a tokamak. Perhaps fusion research, too, might profit from the methods 
developed here. 

5 Retro and prospects 

We have now two systematic approaches to exact solutions of Maxwell's 
equations when material properties vary in space. What was known before 
this article was written? 

The interest in electromagnetically variable media increased dramatically 
with the advent of dielectric waveguides, i.e. in the seventies of the previous 
century. Attempts were made to improve the fibers using dielectrics with a 
graded index of refraction (GRIN), see [TT] for references. In 1975, Kogelik 
noted down a Helmholtz equation with a variable coefficient as in ([2]) which, 
according to his belief, would found a theory of electromagnetic waves in 
GRIN media [lU Sec. 2. 4]. He remarked the similarity of his equation with 
Schrodinger's and rewrote some solutions found in textbooks on quantum 
mechanics for his purpose. For most reseachers this still seems to be the state 
of the art. In 2010, for example, Lipson, Lipson and Lipson struggle to derive 
that Helmholtz equation from Maxwell's equations [151 Sec. 10. 1.2]. They 
encounter a term which bears some similarity with the peculiar second term in 
fISS]) . but cannot handle it and so they discard it pretending an approximation 
called "weak guidance". One can find many examples of proclivity for the 
Helmholtz equation. Maybe for another confirmation of the slow progress, 
see Rauh, Yampolskaya and Yampolskii who reach in 2010 no higher point of 
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view than Kogelik in 1975. They call "Master equation" what others denote 
as Helmholtz or Schrodinger equation. The paper is valuable nevertheless for 
some historical references [19]. 

In their classic monography, Born and Wolf [3[ Sec. 1.6] reproduce some 
deliberations probably first thought by Abeles [1] in 1950. Maxwell's equa- 
tions are written in a cartesian coordinate system and permittivity as well 
as permeability are allowed to depend on one coordinate. Maxwell's system, 
which usually couples all components of the electromagnetic field, is shuffled 
until there is one equation for only one component. This equation has the 
same shape as flMl) although its physical content is different. Unfortately this 
useful equation comes with a second, more complicated equation which has 
to be solved at the same time if the electromagnetic field is to be calculated. 
In summary, the theory put forward by Born and Wolf is practically useless. 

So Born and Wolf haste to a theory which construes the continuum as a 
sequence of small steps. The graded medium is replaced with a pile of thin 
layers. The reflection and the transmission in a single layer are calculated 
from Snell's law and Fresnel's formulae. The results are entered in simple 
matrices such that the reflections and the transmissions in the pile can be 
computed as matrix multiplication. This is the transfer-matrix method. It is, 
with many technical improvements, most popular with practicians. Software 
packages that help to construct the matrices and to execute their multiplica- 
tion can be found and downloaded in the internet. Have a look, for example, 
at "Freesneir'or "RP Coating"! 

In 2010 Turakulov presented a preprint [2^ and in 2011 an article [22] 
wherein vector potentials were used. Ordinary differential equations simi- 
lar to the correct equations fl64l) and (l67I) were found although in a differ- 
ent mathematical context and a useless scalar potential was introduced f2T\ 
which obscures the calculation of the electromagnetic fleld. Turakulov also 
communicated that the equations flMl) and fl671) can actually be solved for 
the example fl85l) . Unfortunately the solutions presented by him were not 
correct. 

The awkwardness of the approaches mentioned so far is caused by the 
arbitrary eliminations. One of the 12 components of the electromagnetic fleld 
and one of 8 Maxwellian equations is selected whereupon lots of unsystematic 
attempts are made to eliminate the other 11 components from the arbitrarily 
selected equation. This is sometimes feasible in the cartesian coordinate 
system. It is cumbersome in cylindrical coordinates and becomes a nightmare 
in spherical coordinates. 

Most physical laws are formulated as partial differential equations be- 
tween between vector flelds and scalars. Thus we have several independent 
variables, denoted by physicists as coordinates, and several dependent ones, 
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denoted by physicists as fields. Given are partial differential equations where 
all these variables are mingled. Wanted are ordinary differential equations 
each with one dependent and one independent variable only. To reach the 
wanted end, we must perform two separations: a separation of the depen- 
dent variables and a separation of the independent ones. These two kinds of 
separation must be kept separate. The fields must be uncoupled without ref- 
erence to special coordinates. This is the idea not comprehended in previous 
work. 

What one must do is this: In a problem with vector fields one must de- 
compose them into their longitudinal and their transverse parts. The former 
can be represented by a scalar potential, the latter by two vector potentials, 
the one being a simple vector potential, the other a vector potential's vec- 
tor potential. Both vector potentials must consist of a predetermined vector 
field times an amplitude which describes the dynamics. Otherwise there is 
no chance to arrive at a differential equation for one scalar quantity only. 

This is the recipe that works in all vector-field theories, e.g. in the theory 
of sound, hydrodynamics, elasticity and electrodynamics [S]. It was applied 
in [7] for a general method solving Maxwell's equations when the material 
constants are constant. Finally it was applied here, namely in (15011511) and 
in 0128111291) . 

The solutions of the equations fl60|) and fl66|) in the Sections l3 . II appear all 
to be new. Seemingly nobody has noticed that the theory of dielectric mirrors 
is based on the Mathieu equation and its generalizations. Calculations up to 
now were done using the transfer-matrix method described above; consider 
the paper by Fink, Winn, Fan, Michel, Joannopoulos and Thomas ^ as an 
anchor of references. Appearently there were not the least precursors for the 
solutions of the spherical problem in Sections 13. 2l ff. nor were there harbingers 
of the two-dimensional theorem of representation in Section HJ especially of 
the solution of the cylindrical problem in Section 14.21 

The broad scope of this article is at the same time its weakness. The 
apphcations were here only indicated, but must be elaborated in detail to be 
useful for the design of optical instruments. For example, in the theory of 
dielectric mirrors, the permittivity must be formulated as a Fourier series 

N 

e= E ene*"'°' (166) 

n=-N 

with positive wavenumber and coefficients e„, n = 0, ±1, ±2, . . . ±A^ which 
model the real medium. Solving the respective equations and (jHS]) costs 
some work, but it can be done as precisely as wanted. The reward are exact 
solutions of Maxwell's equations. These solutions should be better than 
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those obtained using the transfer-matrix method because gradual transitions 
between the layers in the mirror can be taken into account. 
O.Gonzales' hunt for references is gratefully acknowledged. 
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